Decoherence in a classically chaotic quantum system: entropy production and 

quantum— classical correspondence 



O 

o 

(N 

in 



> 
o 
a\ 
o 
\o 
o 



Oh! 

^ : 



Diana Monteoliva * and Juan Pablo Paz ^ 

Departamento de Fisica Juan Jose Giambiagi, FCEyN, UBA, Pabellon 1, Ciudad Universitaria, 1428 Buenos Aires, Argentina 

(February 1, 2008) 

We study the decoherence process for an open quantum system which is classically chaotic 
(a quartic double well with harmonic driving coupled to a sea of harmonic oscillators). We 
carefully analyze the time dependence of the rate of entropy production showing that it has 
two relevant regimes: For short times it is proportional to the diffusion coefficient (fixed by 
the system-environment coupling strength) . For longer times (but before equilibration) it is 
fixed by dynamical properties of the system (and is related to the Lyapunov exponent). The 
nature of the transition time between both regimes is investigated and the issue of quantum 
to classical correspondence is addressed. Finally, the impact of the interaction with the 
environment on coherent tunneling is analyzed. 

PACS numbers: 03.65.Bz 



I. INTRODUCTION 

The process of decoherence has been identified as one 
of the main ingredients to explain the origin of the classi- 
cal world from a fundamentally quantum substrate [1,2]. 
In fact, according to the decoherence paradigm, classical- 
ity is an emergent property imposed upon open systems 
by the interaction with an external environment. In all 
realistic situations this interaction, as it became clear 
in recent years, generates a de facto superselection rule 
that prevents the stable existence of most of the states 
in the Hilbert space of the system. Only a small set of 
states of the system are relatively stable, they are the 
so-called pointer states [3,4]. While pointer states are 
minimally disturbed by the interaction with the environ- 
ment, coherent superpositions of such states are rapidly 
destroyed by decoherence. Thus, this process transforms 
the quantum state of the system into a mixture of pointer 
states. In recent years the study of the physics of deco- 
herence has helped to clarify many interesting features 
of this process. For example: the nature of the decoher- 
ence timescales is now well understood [5,6]; the essen- 
tial features of the process by which the pointer states 
are dynamically selected by the environment are well un- 
derstood [7-9] (see [10] for a recent review). Moreover 
(and most notably) the study of decoherence became ac- 
tive from the experimental point of view where the first 
generation of experiments exploring the fuzzy boundary 
between the quantum and the classical world are already 
starting to produce interesting results [11-13]. 

Over the course of these studies it became clear that 
the decoherence process has very peculiar features for 
quantum systems whose classical analogues are chaotic. 
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In fact, for such systems decoherence seems to be abso- 
lutely essential to restore the validity of the correspon- 
dence principle violated for very short times (the breakup 
time depends logarithmically on the Planck constant) 
[14-17]. The reason for the breakdown of the correspon- 
dence principle for chaotic systems (and for its restora- 
tion due to decoherence) can be understood as follows: 
For these type of systems, quantum evolution continu- 
ously generates a coherent spreading of the wave function 
over large scales both in position and momentum. Thus, 
a classically chaotic Hamiltonian generates a quantum 
evolution which typically produces "Schrodinger cat" 
type states, i.e. starting from a state which is well lo- 
calized both in position and momentum the quantum 
evolution produces a state which is highly delocalized ex- 
hibiting strong interference effects. One may be tempted 
to argue that this effect, while existing, could not be rel- 
evant for large macroscopic systems. But, surprisingly 
or not, even for as large an object as the components of 
the solar system which are chaotic, quantum predictions 
are alarming: On a timescale tft as short as a month for 
Hyperion, one of the moons of Saturn whose chaotic tum- 
bling motion has been analyzed [18], the initial Gaussian 
state of this celestial body would spread over distances 
of the order of the radius of its orbit! [15]. Thus, though 
planetary dynamics appears to be a safe distance away 
from the quantum regime, as a consequence of the chaotic 
character of the evolution, a simple application of the 
Schrodinger equation would tell us that this is not the 
case. In fact, the macroscopic size of a system (be it a 
planet or a cat) is not enough to guarantee its classical- 
ity. Thus, classicality in such system would emerge only 
as a consequence of decoherence, as we will later discuss 
in this paper. 

The reason why a classically chaotic Hamiltonian gen- 
erates highly nonclassical states can be related to the 
fact that chaotic dynamics is characterized by exponen- 
tial divergence of neighboring trajectories. To be able 
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to present this argument, based on the notion of trajec- 
tories, it is better to formulate quantum mechanics in 
phase space (a task that can be accomphshed by using, 
for example, the Wigner distribution [19] to represent 
the quantum state). In fact, if one prepares a quan- 
tum system in a classical state, with a Wigner function 
well localized in phase space and smooth over regions 
with an area which is large compared with the Planck 
constant, it will initially evolve following classical tra- 
jectories in phase space. Therefore, after some charac- 
teristic time the initially smooth wave packet will be- 
come stretched in one (unstable) direction and, due to 
the conservation of volume, squeezed in the other (sta- 
ble) one. This squeezing and stretching is also accompa- 
nied by the folding of classical trajectories which, in the 
fully chaotic regime, tend to fill all the available phase 
space. The combination of these three related effects - 
squeezing, stretching and folding - forces the system to 
explore the quantum regime. Thus, being stretched in 
one direction the wavepacket becomes delocalized and, as 
a consequence of folding, quantum interference between 
the different pieces of the wavepacket, which remain co- 
herent over long distances, develops. The timcscale for 
the correspondence breakdown can be also estimated via 
a simple argument: As the wave packet squeezes expo- 
nentially fast in one direction (momentum, for example), 
ap{t) = ap{0)e~^*, it will correspondingly become coher- 
ent over a distance that can be estimated from Heisen- 
berg's principle as l{t) > ^r^e'^*. When the spread- 
ing is comparable with the scale x where the potential 
is significantly nonlinear, folding will start to apprecia- 
tively affect the evolution of the wave packet and long 
range quantum interference will set in. This time, which 
corresponds approximately to the time when quantum 
and classical expectation values start to differ from each 
other, can therefore be estimated as tn = A~^ln(^m^) 
(note that the quantity x can be as large as the size of 
the orbit of a planet [15,20] for gravitational potential). 

But nothing of this sort is evident in real life: the 
moon of Saturn does not spread over large distances and 
the correspondence principle seems to be valid for macro- 
scopic objects that, like Hyperion, behave according to 
classical laws. So, how can classical mechanics be re- 
covered? Decoherence is a way out of this problem. As 
the system evolves while continuously interacting with 
its environment, it may become classical if the interac- 
tion is such that it continuously destroys the quantum co- 
herence which arc dynamically generated by the chaotic 
evolution. The role of decoherence in recovering the cor- 
respondence principle has been suggested some time ago 
[15] and numerically analyzed more recently [16]. 

But, as originally suggested by Zurek and one of us 
[14], the decoherence process for classically chaotic sys- 
tems has another very important feature, that comes as 
a bonus. The interaction with the environment destroys 
the purity of the system since they become entangled. 
Thus, information initially stored in the system leaks 



into the environment and therefore decoherence is al- 
ways accompanied with entropy production. This is, of 
course, true both for classically regular and classically 
chaotic systems. But what distinguishes chaotic systems 
is the existence of a robust range of parameters for which 
the rate at which the information flows from the system 
into the environment (i.e., the rate of entropy produc- 
tion) becomes entirely independent of the strength of the 
coupling between the system and the environment and 
is dictated by dynamical properties of the system only 
(i.e., by averaged Lyapunov exponents). The reason for 
this can also be understood using a simple (clearly over- 
simplified) argument. The decoherence process destroys 
quantum interference between distant pieces of the wave 
packet of the system and puts a lower bound on the 
small scale structure that the Wigner function can de- 
velop. Thus, the Wigner function can no longer squeeze 
indefinitely but it stops contracting as a consequence of 
the interaction with the environment, which can be typ- 
ically modelled as diffusion [22 24]. As a consequence of 
this, and due to the fact that expansion (or folding) is 
not substantially affected by diffusion, the entropy of the 
system grows at a rate which is essentially fixed by the 
average rate of expansion given by the average Lyapunov 
exponent. In this regime, the entropy production rate 
becomes independent of the diffusion constant (which, 
on the other hand, is responsible for the whole process). 
This result was first conjectured in [14]. More recently, 
numerical evidence supporting the conjecture was pre- 
sented [17,25-28]. The aim of this paper is to present 
solid numerical evidence supporting this result and to 
study other related aspects of decoherence for a particu- 
lar chaotic system. As will become clear later, our studies 
show that the time dependence of the entropy production 
rate has two rather different regimes. First, there is an 
initial transient where the entropy production rate is pro- 
portional to the system-environment coupling strength 
(i.e., in this initial regime the rate behaves in the same 
way as in the regular non chaotic- case). Second, af- 
ter this initial transient the entropy production rate goes 
into a regime where its value is fixed as conjectured in 
[14] by the Lyapunov exponents of the system and is 
independent of the strength of the coupling to the envi- 
ronment. In our work we show that the transition time 
tc between both regimes is linearly dependent on the en- 
tropy of the initial state, and logarithmically dependent 
on the system-environment coupling strength. 

Our aim in this paper is not to give an extensive list of 
all relevant references connected to the study of dissipa- 
tive effects in classically chaotic quantum systems. How- 
ever, we believe it is worth pointing out that our work 
is by no means the first attempt to study these effects 
which, in connection to problems such as the impact of 
noise on localization, the appearence of classical features 
in phase space (like strange attractors, etc) were studied 
in pioneering works several years ago (see for example 
[29]). As general source of reference in this area we rec- 
ommend to interesting books where one can find a rather 
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extended compilation of important works [30,31]. For our 
numerical study we have chosen as nonlinear system the 
driven quartic double well, described in detail in Section 
II. We study the evolution of our system for very simple 
(Gaussian) initial conditions (i.e., initial conditions that, 
from a start, do no exhibit any quantum interference ef- 
fects) so as to focus our attention on the interplay be- 
tween two competing effects: generation and destruction 
of dynamically generated interferences. We also present, 
in Section II, a discussion of the breakdown of correspon- 
dence for chaotic systems. Then, on Section III the way 
in which we model the interaction between our system 
and an environment is presented with some detail. We 
do this by using two complementary approaches based 
on master equations both giving consistent results but 
one of them enable the study of the evolution for long 
dynamical times, which is relevant for discussing the im- 
pact of decoherence on tunneling. Numerical results con- 
cerning the time dependence of the decoherence rate and 
a discussion of their relevance are presented in Section 
IV. Finally, we discuss the impact of the interaction with 
the environment on the tunneling phenomena in Section 
V and we end with some concluding remarks on Section 
VI. 



II. EVOLUTION OF THE QUANTUM ISOLATED 
SYSTEM 

In this Section we study the dynamical behaviour of 
an isolated quantum system with Hamiltonian 



of the driving amplitude, chaotic layers are so thin that 
only those near the unstable fixed point of the integrable 
system can be seen in the stroboscopic phase space por- 
trait. As the value of s is increased these chaotic regions 
gain in relevance while resonances other than the two first 
ones near the stable fixed points become in turn smaller 
and are almost invisible to the eye. For large values of 
the driving amplitude s the different chaotic layers merge 
into what is called the chaotic sea; the stable islands and 
two first resonant islands are much reduced (see Figure 
la). In what follows we will choose parameters so that 
either the majority of the phase space is chaotic, as in 
Figure lb, where the stable islands around the stable 
fixed points of the integrable system have shrunken to 
invisibility and only the two first resonance islands are 
seen, or the regular regions coexist with the chaotic sea 
as shown in Figure la. 




Ho{x,p,t) = bx'^ + — 1- sxcos{u)t), (1) 

which corresponds to a harmonically driven quartic dou- 
ble well. This non-linear system has been extensively 
studied [32,33] in the literature. For s = the system 
is integrable and exhibits the generic phase space of a 
bistable system. It has two stable fixed points &t x = 
±y^326a), p = (with associated energy E = — 166^a) 
and one unstable fixed point at x = 0, p = (with 
E = 0). The presence of the driving term introduces an 
infinite number of primary resonance zones in the phase 
space of the system. The tori near the separatrix of the 
non-driven system {E = 0) and KAM invariant tori in 
the region between resonance zones are progressively de- 
stroyed as the amplitude s of the driving increases. The 
phase space of the system starts to look chaotic, present- 
ing in general a mixed nature as can be seen (for two 
sets of parameters) in Figure 1. The regions around the 
two stable fixed points and the chains of resonant islands 
form regular regions where the dynamics of the system 
is nearly integrable. The chaotic layers between regular 
regions develop from the homoclinic tangles, the intrin- 
cated interweaving of the stable and unstable manifolds 
originated at the hyperbolic fixed points between the res- 
onant islands or the one at the origin. For small values 
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FIG. 1. Stroboscopic phase space of the driven double well 
with parameters (a) m = 1, 6 = 10, a = 1/32, s = 1, ui = 5.35 
and (b) m = 1, 6 = 10, a = 1/32, s = W, ui = 6.07. The 
borders of the dark elipses represent contours of minimun un- 
certainty Gaussian initial conditions at 1/20 of its peak value, 
for xq — —3.7, po = 0, Ux ~ .05, = 1 on the leftmost island 
and xo = 1, po = 0, = .05, C7p = .05 on the chaotic sea. 

As we mentioned above, our aim is to follow the quan- 
tum evolution of initial states whose Wigner functions 
are fully contained either in a regular island or in the 
chaotic SGcl^ clS illustrated in Figures la and lb. This re- 
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quircmcnt forces us to take values of h which arc small 
enough. Indeed, for the parameter set corresponding to 
Figure la, the leftmost regular island has an area which 
is approximately given by ^ ~ Qn. If we need states 
well localized within each region (the regular islands or 
the chaotic sea) we need to compare this area with the 
one covered by the dark ellipses appearing in Figure 1, 
which are given by Ae = 2'!T\n{20)ax(Jp = 7rln(20)?i (i.e., 
we consider initial states which are pure, minimal un- 
certainty, coherent states). Thus, a reasonable value for 
U (the one we choose) is h = 0.1 for the island to be 
about 20 times larger than the extent of the initial state 
(for the parameters corresponding to Figure lb a larger 
value of h could also be chosen for our purposes). + 
We numerically solved the Schrodinger equation for this 
system evolving the above initial states using two differ- 
ent methods. First, we integrated this equation with a 
step by step algorithm using a high resolution spectral 
method [34]. Second, we solved the same equation ap- 
plying a numerical technique based on the use of Floquet 
states [35]. This method consists of numerically find- 
ing eigenstates of the unitary evolution operator for one 
period (Floquet states) which form a complete basis of 
the Hilbert space of the system. By expanding the ini- 
tial state on this basis, the solution of the Schrodinger 
equation becomes trivial [36]. Thus, all the difficultly of 
the method is hidden in finding the Floquet states (the 
same applies when solving the Schrodinger equation for 
a time-independent Hamiltonian). We succesfully used 
this method for some parameter sets but did not apply it 
for all cases of interest. In fact, the main difficulty of the 
method is the number of Floquet states required to accu- 
rately expand the initial states shown in Figure 1, which 
rapidly grow as h becomes small (for example, for the 
parameters of Figure lb, with h = 0.1, one needs more 
than 200 Floquet states to evolve the wave packet cen- 
tered in the chaotic sea). In what follows we will discuss 
the typical results one finds when solving the Schrodinger 
equation with the above initial states. 

A. Initial states in a regular island 

The time series of the expectation value of the po- 
sition X of the particle both for classical and quantum 
evolutions, for the initial Gaussian state located within 
the leftmost regular island as shown on Figure la can 
be seen in Figure 2. Notice that the state is cen- 
tered at Xo = — 3.7,po = and that its variances are 
= .05,f7p = 1. Classical and quantum expectation 
values remain identical within the numerical errors of our 
calculations, in agreement with previous results [33,37]. 
Similar results were found for higher order moments of x 
and p and for different initial conditions embedded either 
in the stable or the resonant islands and for different sets 
of parameters of the system. These results can be under- 
stood as follows: When the initial Gaussian function is 



located within a regular island, its Floquet decomposition 
is characterized by a small number of states, mostly local- 
ized on the regular regions. Thus, only very few frequen- 
cies enter into play in the dynamics and the time series of 
the expectation value of any observable appears as quasi 
periodic and regular. Moreover, as the Floquet states 
entering in the decomposition of the initial state are lo- 
calized on integrablc tori, EBK quantization is possible 
and both classical and quantum expectation values (and 
distribution functions) look identical for long time scales. 
Breakdown of correspondence is expected on time scales 
inversely proportional to some power of h [38], which is 
sufficiently slow to cause no difficulties with the classical 
limit of quantum theory. 
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FIG. 2. Time scries of the classical and quantum expecta- 
tion value of X for the states drawn on Fig. la. 

Nevertheless, as localized Floquet states come from 
the superposition of almost degenerate symmetry-related 
pairs (doublets), each belonging to a different par- 
ity class, the quantum particle will eventually tunnel 
through the chaotic sea into the related regular island 
[33,39,40] while for the classical particle this feature is 
forbidden. Tunneling times are very long for the param- 
eters we chose (typically they are of the order of 50-100 
times the driving period). The accumulation of error 
makes step by step integration of the Schrodinger equa- 
tion troublesome to study this tunneling regime. How- 
ever, a solution based on Floquet states is more con- 
venient since it allows us to compute the state at any 
time provided we know it within the first period r, with 
T = 27r/a; [40-42]. We applied this method for several 
sets of parameters chosen in such a way that the num- 
ber of Floquet states needed to faithfully represent the 
Hilbert space doesn't become too large (still enabling us 
to keep states well localized within the regular island). 
One of the parameter sets we used for this purpose is 
m = 1, 6 = 10, a = 1/32, s = 4, a; = 5.35 and h = 1. 
The stroboscopic phase space looks like something in be- 
tween Figures la and lb but the value of h is ten times 
larger than the one used in the above Figure. For this 
parameter set, Figure 3 shows the evolution of the ex- 
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pectation value of the position for an initial state located 
in the leftmost regular island, centered at xq = —3.52, 
Po ~ (with ax = 0.25, (jp — 2). This initial state was 
expanded by 20 Floquet states. In this case we clearly see 
the system tunnel from one regular island to the other one 
(deviation from classical behavior is observed for times 
well below the tunneling time, which in this case of the 
order of t = 56t). 




40 80 t 

FIG. 3. Time series of the quantum (full line) and classical 
(dotted line) expectation value of x for an initial state located 
within a regular island (see text). 

The evolution of quantum and classical distribution 
functions is shown in Figure 4 for two different times: a 
time half way the tunneling {t = 28t) where the wave 
packet is completely delocalized, and the time for which 
the state has tunneled to its pair related regular island 
on the right {t — 56r). Notice that though the state is 
initially peaked in the leftmost regular island, there is a 
nonnegligible probability of finding it on the chaotic sea, 
due to the relatively high value of h. For this reason, 
the tail of the initial Gaussian distribution spreads over 
the whole chaotic sea as time goes on. For the classi- 
cal case, this effect is not very important since the state 
stays within the regular island. On the contrary, for the 
quantum evolution, one clearly observes the wave packet 
tunneling through the chaotic sea into the pair related 
regular island on the right. 
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FIG. 4. The classical (left) and quantum (right) evolution 
of the phase space distribution function for the tunneling case. 
Parameters are that of Fig. 3. 



B. Chaotic states 

For Gaussian states initially located within the chaotic 
sea, as the one corresponding to the ellipse in Figure la, 
the behaviour is quite different from that of the regular 
initial state. In fact. Figure 2 shows the time series of 
the quantum and classical expectation values of the po- 
sition for this initial state. We see that after a time scale 
much shorter than that of the regular case (our results 
are consistent with a logarithmic dependence of this time 
scale on h) , discrepancies between classical and quantum 
predictions are evident. For the various different initial 
conditions and parameters of the system tested, break- 
down of correspondence occurs at very early times, and 
in all cases it is comparable with the dynamical scale 
T. Apart from the noticeable deviation between quan- 
tum and classical predictions one clearly sees that the 
expectation values have a much more irregular behav- 
ior than in the integrable case. This is expected and, in 
the quantum regime, can be explained as due to the fact 
that the number of Floquet states required to expand 
the initial state is, contrary to the regular case, rather 
large. Thus, as Floquet states are mostly extended over 
the whole chaotic sea (all the more in the semiclassical 
limit) the evolution of the initial Gaussian state will con- 
tain many more frequencies (a few hundred, typically) 
and as a consequence of this the time evolution of any 
expectation value looks rather irregular. The deviation 
between classical and quantum expectation values is just 
a consequence of the fact that the quantum state, de- 
scribed in the phase space by the Wigner function, be- 
comes more and more different from the corresponding 
classical distribution function [16]. This situation is illus- 
trated in Figure 5 where we show the Wigner distribution 
function for the above initial Gaussian state (localized in 
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the chaotic sea), compared with the corresponding clas- 
sical distribution function for four different times. Very 
quickly - even before classical and quantum expectation 
values deviate from each other -the quantum distribu- 
tion function becomes affected by small scale interference 
effects, and develops into a highly non-classical distribu- 
tion function. This results are also in agreement with 
previous ones [33,37]. 
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FIG. 5. Classical (left) and quantum (right) distribution 
functions for four different times. Already for early times 
interference effects are observed, but after t ~ 2r they sprout 
over the whole available phase space and go down to very 
small scales, as shown by the insets of area h, which blow-up 
the tiny rectangles drawn on the Wigner functions for t — 2t 
and t = 9r. 
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In this equation, the first term on the right hand side is 
the Poisson bracket that generates purely classical evolu- 
tion. The second term, containing higher odd derivatives, 
is responsible for the quantum corrections [19,37]. The 
fate of the Wigner function can be understood qualita- 
tively by using the following simple argument [14]: If one 
starts from an initial state which is smooth (as the ones 
shown in Figure I), the higher derivatives terms are negli- 
gible and the Poisson bracket term dominates. Therefore, 
the Wigner function will tend to follow initially the Li- 
ouville flow, i.e., it will evolve following the classical tra- 
jectories. As any point in the chaotic sea is hyperbolic, 
the initially regular patch will stretch along the unstable 
manifold and squeeze in the other - stable - direction. 
As the Wigner function squeezes, its derivative will tend 
to grow and, consequently, quantum corrections will tend 
to become more and more important. The time scale for 
which quantum corrections must be taken into account 
can be estimated as tn = \n{xcrq{0)/h) [14] where x is 
the scale where the nonlinearities of the potential come 
into play and is usually defined as x = iOxV/d^Vy/"^ 
{<Jq{0) is the initial spread of the distribution function 
in the stable manifold and A is the Lyapunov exponent 
of the system). Moreover, as the motion is bounded, the 
Wigner function will eventually tend to fold and so differ- 
ent pieces will coherently interfere making the distribu- 
tion function develop small scale structure as seen in Fig- 
ure 5. The scale at which structure will tend to develop is 
typically sub-Planckian: thus, as we will have Sp — h/L 
and 6x = Ti/P (where L and P are the size of the system 
in position and momentum), the region over which the 
Wigner function will tend to oscillate has an area approx- 
imately given hy Ara / LP = hx (h/LP) < h [43]. In 
our simulations we verified that, after a time consistent 
with tn ~ 0.5r, for the set of parameters here considered, 
the quantum phase-space distribution does not resemble 
the classical one and the Wigner function oscillates wildly 
on tiny scales of size 6p « .01 and Sx « .005 (being the 
product of the order of ^ = 6p5x w 5 x 10^'"' ^ fi). The 
insets on Figure 5 show the development of small scale 
structure on the Wigner functions for relevant times. 
There, a detail of the distribution function over a region 
of area h is shown. The existence of this sub-Planckian 
structure in the Wigner distribution has been overlooked 
in the literature and its relevance has only been noticed 
recently [43]. 



For our purposes it is very useful to represent the tem- 
poral evolution of quantum states in phase space. The 
time dependence of the Wigner function (shown in Fig- 
ure 5), is goverened, in general, by an equation entirely 
equivalent to Schrodinger equation, which can be written 
as: 



III. COUPLING TO THE ENVIRONMENT: 
MASTER EQUATIONS 

Here we will describe the way in which we analyze the 
evolution of our system when it is coupled to an cnviron- 



6 



ment. The role of such environment will be played in our 
case by an infinite number of oscillators coupled to the 
system via an interaction term in the Hamiltonian, which 
is assumed to be bilinear both in the coordinates of the 
system and the oscillators. As we are interested in follow- 
ing the; c;volution of the system solely, we will compute the 
reduced density matrix p obtained from the full density 
matrix of the Universe (system+environment) by tak- 
ing the partial trace over the environment p = TrE{pu)i 
where pu is the full density matrix of the Universe. This 
reduced density matrix will obey a master equation which 
can be obtained from the full von-Neumann equation by 
tracing out the environment. There is a vast literature 
dealing with the properties of the master equation ob- 
tained under a variety of assumptions. Our purpose here 
is not to present a review of these results but to sketch 
the basic method used to obtain the master equations 
we will use for our studies below (see [10] for a more ex- 
tensive review on the derivation and properties of master 
equations in studies of decohcrcncc) . 

The general approach to obtain the master equation 
is the following: We model the environment by a set of 
harmonic oscillators with mass and natural frequency 
LOi, in thermal equilibrium at temperature T. The Hamil- 
tonian is Hn = X^jWj {blbi + 1/2), where bi and &J are 
the annihilation and creation operators, respectively, for 
a boson mode of frequency uji. We further assume that 
the interaction Hamiltonian between the system and the 
environment is i?/ = x J2i{9ibi + 9ibl), where gi are 
coupling constants. 

The civolution of the reduced density matrix of the par- 
ticle can be obtained by taking the partial trace over the 
environment of the exact von Neumann equation for the 
full Hamiltonian, that reads i h p'u = [ H(t) , pu] . This can 
be straightforwardly done under a number of standard 
assumptions (see [10] for more details): First, one as- 
sumes that the system and the environment are initially 
uncorrelated and that the reservoir is in an initial state of 
thermal equilibrium at some temperature T. Second, one 
assumes that the system and the environment are very 
weakly coupled. So, solving the von Neumann equation 
perturbatively in the interaction picture (up to second 
order in perturbation theory) a master equation for the 
reduced density operator is obtained, which fully deter- 
mines the quantum dynamics of our system. It reads: 

p{t) = - I dt' gs{t-t')[x{t),[x{t'),p{t)]] + 
Jo 

+ 1 [ dt' QAit-t') [X{t),{x{t'),p{t) }] , (3) 

where the above kernels are 

9s{t -t')= |ffi|2(l + 2n{uJi)) cos{uJi{t - t')) , 

gA{t-t')= J2\gi\^sm{u)i{t-t')) , (4) 



with n{uj) = ^i3ht_i (and /? = l/fc^T). It is worth stress- 
ing that the above master equation is derived without 
appealing to the usual Markovian approximation (see 
[10,22]). However, in our studies we restrict ourselves to 
the Markovian regime where the kernels (4) are local in 
time. We solved the above mastcir equation for two differ- 
ent regimes and environmental couplings. First, we con- 
sidered the widely used ohmic-high temperature regime 
of the Brownian motion model where a simple master 
equation can be written and solved. Second, we obtained 
a master equation for the evolution of a coarse grained 
density matrix obtained from eq. (3) by averaging over 
one driving period. This can be naturally done by using 
the Floquet representation and is useful to study some 
properties of the open system in the long time regime. 
We will now describe the two methods. 



A. Ohmic— high temperature environment 

The simplest special case of equation (3) follows when 
assuming the high temperature limit of an Ohmic envi- 
ronment. Defining an spectral density for the environ- 
ment as J{ui) = limAw^o Ea><w*<w4-Au; IffiP, assum- 
ing that the spectrum is ohmic, i.e., that J{uj) = ^muj/h, 
and going to the high temperature limit, gs{t — t') = 
D5{t-t') and gA{t-t') = 2jS{t-f) with D = 2m'jkBT, 
the following well known equation for the Wigner func- 
tion of the system is obtained: 

w = {Ho, w},, + x: 2i(2;+i), ^^'"^^^^^r^^^^ 

+ 2^dp {pW) + Dd^pW. (5) 

The physical effects included in this equation are well 

understood: The first term on the right hand side is the 
Poisson bracket generating the classical evolution for the 
Wigner distribution function W; the terms in h add the 
quantum corrections. The environmental effects are con- 
tained in the last two terms generating, dissipation and 
diffusion, respectively. In the high-T limit here consid- 
ered dissipation can be ignored and only the diffusive 
contributions need be kept (doing this we are taking the 
7^0 limit while keeping D constant). The diffusion 
constant was chosen small enough so that energy is nearly 
conserved over the time scales of interest. Equation (5) 
was integrated by means of a high resolution spectral al- 
gorithm [34]. Results were stable against changes both 
in the resolution of phase space and the time accuracy 
required on each integration step. They will be presented 
on Section 4. 
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B. Coarse grained master equation in the Floquet 

basis: 

It is interesting to write equation (3) in the Floquet 
basis \4'n{t)} = exp"*'^''* \(j)^{t)}, where the \(p^i{t)) are r- 
periodic, taking advantage of the r-symmetry of our sys- 
tem Hq. There are in the hteraturc several approaches 
of the kind [39-42]. Ours is similar to that of Kohlcr 
et. al. [44], though in that reference the authors have 
restricted the use of the equation to the study of chaotic 
tunneling near a quasienergy crossing. When one writes 
the master equation in the Floquet basis and takes the 
Markovian high temperature limit, the resulting equation 
has T-periodic coefficients. Therefore, one can derive a 
temporal coarse grained equation by taking the average 
of such equation over one oscillation period. The result- 
ing equation for the average density matrix in one period 
a turns out to be 

0-^,1, = '^M^^v,a,l30'a,l3, (6) 

where the coefficients M^^^^a,/} are defined as 

+ SaA{^^))0,. - 2 {{X))0,.} (7) 

where the notation = i 

is used to denote time averages of matrix elements of op- 
erators in the Floquet basis. As this is an equation with 
constant coefficients, once these coefiicients are numeri- 
cally calculated, the solution is formally obtained for all 
times as = Ea,/3 i^""') ^,.,c,i3^-A0)- The difii- 

culty in solving this equation resides on the large number 
n of Floquet states one typically needs to accurately rep- 
resent a state located in the chaotic sea (as we mentioned 
above, a semiclassical argument to estimate the number 
of such states is given by the ratio A/(2TTh), where A 
is the area of the chaotic sea, which might quickly be- 
come a very large number). As M has dimension n*, 
numerical limitations spring up at this point. However, 
for some parameter values it was still possible to manage 
the numerical problem. The results will be presented on 
Section V. 



IV. RESULTS 

Before going into a detailed description of our model 
for the open system it is instructive to analyze equation 
(5) so as to have an intuitive idea of what is going on. 
When the diffusion term is absent that means, states 
evolve simply according to Schrodinger equation- for a 
smooth initial state the dominant term is the Poisson 
bracket. As wc already discussed in Section 2, the Wigner 
function initially evolves following nonlinear classical tra- 
jectories, looses its Gaussian shape and develops tendrils 



while folding. If the initial state is located in the chaotic 
sea, this happens exponentially fast. Due to the combi- 
nation of squeezing-stretching and folding the gradients 
increase and quantum corrections in equation (5) become 
important, so discrepancies between quantum and clas- 
sical predictions begin to be relevant. Also, quantum 
interferences among different pieces of the wave packet 
develop and generate oscillations in the Wigner function. 

The effect of the dccoherence producing term in equa- 
tion (5) can be understood as being responsible of two 
interrelated effects: On the one hand, the diffusion term 
tends to wash out the oscillations in the Wigner function 
suppressing quantum interferences. Thus, for this system 
decoherence is the dynamical suppression of the interfer- 
ence fringes that are dynamically produced by nonlinear- 
ities. The time-scale characterizing the disappearance of 
the fringes can be estimated easily using previous results 
[6]: Fringes with a characteristic wave vector (along the 
p-axis of phase space) kp decay exponentially with a rate 
given by Td = Dkp. Noting that a wave packet spread 
over a distance Ax with two coherently interfering pieces 
generate fringes with kp = Ax/U one concludes that the 
decoherence rate is F/j = DA^/h . This rate depends 
linearly on the diffusion constant. When the Wigner 
function is coherently spread over the whole available 
phase space one expects fringes with wavelength of the 
order of 1 /L where L is the size of the system (these are 
responsible for the existence of sub-Planck structure in 
the Wigner function, as mentioned in Section 2). The 
rate at which these fringes disappear due to decoherence 
is then To = DL'^/h^. Thus, we can derive a condition 
that should be satisfied in order for these fringes to be ef- 
ficiently destroyed by decoherence: the destruction of the 
fringes, that takes place at a rate F^, should be faster 
than their regeneration, that takes place at a rate f2/ 
fixed by the system's dynamics and related both to the 
Lyapunov time and the folding rate. Thus, if Clj -^Tu 
the small scale fringes will be efficiently washed out and 
the Wigner function would remain essentially positive. It 
is worth mentioning that the destruction of fringes would 
generate entropy at a rate which, provided the condition 
flf <C Tjj is satisfied, should be independent of the dif- 
fusion constant, and should be fixed by f2/. Thus, one 
could argue that every time the Wigner function stretches 
and folds, becoming an approximate coherent superposi- 
tion of two approximately orthogonal states, the destruc- 
tion of the corresponding fringes should generate about 
one bit of entropy. If the timescale for the fringe disap- 
pearance (To) is much smaller than the one for producing 
the above superposition (fixed by fi/) then the entropy 
production rate would be simply equal to f2/ (ideally, 
one would expect one bit of entropy created after a time 
!/%)• 

However, the disappearance of the interference fringes 
is not the only effect produced by decoherence. There 

is a second related consequence of this process (which 
is also present in the classical case): While interference 
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fringes arc being washed out by decoherence, the dif- 
fusion term also tends to spread the regions where the 
Wigner function is possitive, contributing in this way to 
the entropy growth. But, as discussed in [14,28], the rate 
of entropy production distinguishes regular and chaotic 
cases. For regular states, decoherence should produce en- 
tropy at a rate which depends on the diffusion constant 
D. However, for chaotic states the rate should become 
independent of D and should be fixed by the Lyapunov 
exponent. The origin of this D-independent phase can be 
understood using a simple minded argument (presented 
first in [14] and later discussed in a more elaborated way 
in [28]): Chaotic dynamics tends to contract the Wigner 
function along some directions in phase space competing 
against diffusion. These two effects balance each other 
giving rise to a critical width below which the Wigner 
function cannot contract. This local width should be 
approximately cr^ = 2D/\ (being A the local Lyapunov 
exponent). Once the critical size has been reached, the 
contraction stops along the stable direction while the ex- 
pansion continues along the unstable one. Therefore, 
in this regime the area covered by the Wigner function 
grows exponentially in time and, as a consequence, en- 
tropy grows linearly with a rate fixed by the Lyapunov 
exponent. Moreover, the appearance of a lower bound 
for the squeezing of the Wigner function make quantum 
corrections unimportant and so the Wigner function will 
evolve classically (following Liouville flow plus diffusive 
effects), and the correspondence will be reestablished. 

In this Section we present solid numerical evidence sup- 
porting the existence of this D-indepcndent phase. For 
simplicity, instead of looking at the von Neumann en- 
tropy TivN = — Tr {pr log Pr) we examine the linear en- 
tropy, defined as Ti = — log(Tr (p^)), which is a good 
measure of the degree of mixing of the system and sets a 
lower bound on Hvn (i-e, one can show that Hv w > W) . 
The above argument concerning the role of the critical 
width cjc may appear as too simple but captures the es- 
sential aspects of the dynamical process. Wo can present 
a more elaborate argument using the master equation 
(5) to show that the rate of linear entropy production 
can always be written as: 

Ti = 2D{{dpWf)/{W^), (8) 

where the bracket denotes an integral over phase space. 
The right hand side of this equation is proportional to the 
mean square wave-number computed with the square of 
the Fourier transform of the Wigner function. This im- 
plies that the entropy production rate is closely related to 
the phase space structure present in the Wigner distribu- 
tion. Thus, the D-independent phase begins at the time 
when the mean square wave-length along the momentum 
axis scales with diffusion as (as CTc does). This behav- 
ior cancels the diffusion dependence of Ti. which becomes 
entirely determined by the dynamics. 

Appart from analyzing the D independent phase of en- 
tropy production we analyze the nature of the the tran- 
sition between the diffusion dominated to the chaotic 



regime. This time tc can also be estimated along the 
lines of the previous argument: The time for which the 
spread of the Wigner function approches the critical one 
is tc « log(crp(0)/crc). According to this estimate 
tc should depend logarithmically on the diffusion con- 
stant and on the initial spread of the Wigner function 
(for Gaussian initial states the spread depends exponen- 
tially on the initial entropy, therefore tc should vary as 
a linear function of the initial entropy). Our mimerical 
work is devoted to testing these intuitive ideas. 

In the following subsections we present our results. 
First, we will show (for completeness of our presenta- 
tion) how decoherence restores classicality washing out 
interference fringes. Then, we focus on our main goal: 
the study of the entropy production rate of the system 
as a function of time. 



A. Correspondence principle restored: 
disappearance of interference fringes 

In Section 2B we showed how classical and quantum 
expectation values of the position observable become dif- 
ferent from each other after a relatively short time when 
the initial state is located within the chaotic sea. In Fig- 
ure 6 we show how this result is affected by decoherence. 
There one observes the time dependence of the expec- 
tation value of position obtained by solving the master 
equation (5) (i.e., considering the decoherence effect) as 
compared with the corresponding classical time series. It 
is clear that, in accordance with results previously ob- 
tained in [16], the correspondence principle has been re- 
stored. 

5 I • . . . 1 

{x{t)) 




I . . . . 1 

2 4 6 8 i 10 

FIG. 6. Time series of the classical and quantum expecta- 
tion value of X for the states drawn on Figure la, for the open 
system {D = 0.01) 

With respect to the behaviour of the Wigner function 
(shown for the case of pure Schrodinger evolution in Fig- 
ure 3), the impact of decoherence on this distribution is 
seen in Figure 7 where we compare the decohered Wigner 
function with the classical distribution at three different 
times. It is clearly seen that decoherence in this case 
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is strong enough to induce classicality, which is reflected 
at the level of the Wigner function and, consequently, at 
the level of any expectation value. It is worth mentioning 
that the condition for classicality discussed above [14,15], 
i.e. tTcX ^ is satisfied in our case: For the parame- 
ters we are using here we have that CTc ~ 0.2, x ~ 4 
and therefore we are in a (not so highly) classical regime 
(remember that we are using h = 0.1). The small scale 
structure developed by the distribution fuction on the 
isolated evolution (Figure 5) is stopped when it reaches 
the lower bound CTc imposed by the environment. The in- 
sets on Figure 7 shows the portion of area A = of the 
decohered Wigner function marked by the tiny rectangle 
for t = 9r. The sub-Planckian structure observed for the 
isolated evolution is now absent: the interference fringes 
have a typical size which is above Uc = 0.2. 




-1 -2 □ 2 1 -1 -2 □ 2 1 




FIG. 7. Classical (left) and quantum (right) distribution 
functions for the same initial condition as on Fig. 5, at t = 2r 
and t — 9t, when the system is opened to the action of the 
environment {D — 0.01). 

B. Entropy Production 

Here we study the time dependence of the entropy pro- 
duction rate. In Figure 8 we plot the time dependence 
of log(^) both for the initial condition centered in the 
regular island and for the one centered in the chaotic 
sea (the initial states are the ones corresponding to those 
shown in Figure la). Figure 8 illustrates one of the main 
points we want to establish in this paper. First of all 
one notices a drastic difference between the behaviour of 
the entropy production rate for regular or chaotic cases 
. For regular initial conditions, the entropy is always 
produced at a rate which is linearly dependent on the 
diffusion coefhcient D as it is clearly seen in the plot at 
the top of Figure 8. On the other hand, for chaotic ini- 
tial conditions the behavior is completely different. For 
early times the rate depends linearly on D but this initial 



regime is rapidly followed by one where the entropy pro- 
duction rate is independent of the value of the diffusion 
coefficient. The existence of the initial D dependent tran- 
sient for the chaotic case comes as no surprise. Indeed, 
this is what is expected if the entropy is coming from: 
(i) the destruction of interference fringes (which are ini- 
tially generated at a relatively slow rate), (ii) the slow 
increase in the area covered by the possitive part of the 
Wigner function. The oscillations evident in both regu- 
lar and chaotic evolution of the rate have the frecuency 
of the driving force and are related both to changes in 
orientation of the fringes (decoherence is more effective 
when fringes are aligned along momenta) and, more im- 
portantly, to the change in spread of the Wigner function 
in the momentum direction induced by the dynamics. 




5 10 




5 10 15 



time 

FIG. 8. Entropy production rate (in logarithmic scale) vs. 
time (in units of the driving period). The bold curve is the 
(time dependent) Lyapunov exponent. The linear dependence 
of the rate on D appears in the graph at the top (regular ini- 
tial state) and during the initial transient in the lower plot. 
In this case (initial state in the chaotic sea) the rate becomes 
independent on diffusion and is equal to the Lyapunov expo- 
nent (if D is not too small, see text). 

Contrasting with this initial behaviour described 
above, initial conditions on the chaotic sea undergo a sec- 
ond, very different regime, where the entropy production 
rate is independent of the value of the diffusion coefficient 
D. Moreover, the numerical value of the rate oscillates 
around the average local Lyapunov exponent of the sys- 
tem. This is the bold curve shown in Figure 8 that is 
simply computed as the average Lyapunov over an en- 
semble of trajectories weighted by the initial distribution. 
For each trajectory the time dependent local Lyapunov 
exponent was calculated using the method proposed in 
[45]. This result, which is robust under changes on ini- 
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tial conditions and on the parameters characterizing the 
dynamics, confirms the conjecture first presented in [14]. 

It is interesting to remark that while the simple pic- 
ture presented in [14] is in good qualitative agreement 
with our results, the arguments presented in that pa- 
per are too simple to include some important effects we 
found here. In particular the oscillatory nature of the 
rate was completely overlooked in [14]. However, having 
said this, it is still possible to test some important re- 
sults obtained in [14] for the transition time tc between 
both regimes. First, we analyzed the dependence of the 
transition time on the diffusion coefficient D. Due to the 
oscillatory nature of the rate, there is some ambiguity 
in the definition of tc- Here, we defined it as the time 
for which the rate reaches some value after the initial 
transient. As the rate goes through a jump of two or- 
ders of magnitude when changing from one regime to the 
other, this definition is a reasonable one. Thus, we found 
a logarithmic dependence of the transition time on the 
diffusion coefficient, as can be seen on Figure 9. Second, 
we investigated the behaviour of the rate as a function 
of the entropy of the initial state. Our definition of tc is 
the same as before. Parameters of the system for these 
studies where those that would allow states with initial 
entropy up to H(0) = 4 be easily located within the 
chaotic sea. We obtained thus the results shown on Fig- 
ure 9, where a linear dependence of the transition time 
tc on H{0) is clearly seen. Both results confirm the naive 
expectation concerning the nature of tc that we discussed 
above. 

It is remarkable that for long times the entropy produc- 
tion rate is indeed fixed just by the dynamics, becoming 
independent of D (after all, the entropy production is 
itself a consccucncc of the coupling to the environment 
but the value of the rate becomes independent of it!). 
The results presented in Figures 8 and 9 were shown to 
be robust under changes of initial conditions and other 
parameters characterizing the classical dynamics. 
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FIG. 9. The transition time between the diffusion domi- 
nated regime and the one where the entropy production rate 
is set by tlic Lyapunov exponent is shown to depend linearly 
on entropy (top) and logarithmically on the diffusion constant 
(bottom) . Numerical results were obtained using the param- 
eters: B = W,C = 0.5, E = W,Lj = 6.16, D = IQ-^ (top), 
7^(0) = (bottom) 

There are two limitations for the above results to be 
obtained. On the one hand the diffusion constant cannot 
be too strong: In that case the system heats up too fast 
and the entropy saturates, making the numerical simula- 
tions unreliable. On the other hand, diffusion cannot be 
too small either: If that is the case decoherence may be- 
come too weak and the interference fringes could persist 
over many oscillations; the minimal value of D required 
for efficient decoherence could be estimated as described 
above: If the Wigncr function is coherently spread over a 
region of size A^, ~ L, we would need a diffusion constant 
larger than Dmin ~ fi' jL"^ ^ 10~^ for the environment to 
be able to wash out the smallest fringes in one driving pe- 
riod (this is nothing but the above condition flf <C Tjj). 
In fact. Figure 8 shows that when D = 10~^ the entropy 
production rate is one order of magnitude smaller than 
the one corresponding to D = 10~^. Thus, for values of 
D that are too small,the condition for classicality is not 
satisfied and the D independent phase of the evolution 
is never attained: the Wigner function always retain a 
significant negative part and decoherence is not effective. 



V. DECOHERENCE AND THE SUPPRESSION 
OF TUNNELLING 

For initial states localized in the regular islands corre- 
spondence is broken for long times when tunnelling be- 
comes effective, as pointed in Section 2A. Here, we in- 
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vestigate the influence of decoherence on this process by 
using our coarsed-grained master equation (6) . In doing 
this we should be careful! to choose a number of Flo- 
quet states which is large enough. Thus, decoherence 
couples these states and the expected quasi-equilibrium 
state one gets from the master equation should be ap- 
proximately diagonal in the Floquet basis. One expects 
that using n Floquet states to expand our Hilbert space, 
the master equation would tend to mix the state in such 
a way that the entropy would grow up to a level where all 
states become occupied with equal probability. At this 
point the numerical simulation done in this way becomes 
clearly unreliable. We solved the master equation com- 
puting the Von Neumann entropy Hvn = —Tr {a log a) 
checking that its value is kept well below saturation. For 
the same parameters we described in Section 2A (and us- 
ing Floquet 40 states instead) we were able to accurately 
study the tunnelling process from an initial state local- 
ized in the leftmost regular island. We analized the evo- 
lution of the expectation value of position (x) — Tr[xa] 
which is shown in Figure 11. Suppression of tunneling 
is clearly observed. Notice that the asymptotic value Xa 
of (x), though small, is not zero. This value can be es- 
timated as Xa = limj^oo /^^ p^^j^ (ia;dpa;W^(a;,p, where 
is the region of the phase space within the leftmost 
regular island where the state is initially located. The 
numerical result shown in the above Figure is consistent 
with this estimate, meaning that the final state a has a 
significant part which remains trapped in the left of the 
well. The evolution of the Wigncr distribution function, 
shown on Figure 10 for t = 28t and t — 56r, illustrates 
the described behaviour of the state of the system. It is 
also noticeable that the state is trapped in the leftmost 
island as a result of the interaction with the environment. 
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FIG. 10. Classical (left) and quantum (right) evolution of 
the Gaussian state initially localized at the leftmost regular 
island (see Section 2A) when the system is opened to the 
action of the environment environment {D = 0.01), for the 
tunneling case. 
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FIG. 11. Time series of the quantum expectation value of 
X for the state of Fig. 3 

It is interesting to notice that this picture would dras- 
tically change if the coupling to the environment is 
switched on at a time > (i.e., if we let the sys- 
tem to evolve freely, with its own Hamiltonian, before 
coupling it to the environment). This is very easy to do 
using the above master equation. The results obtained 
in this way are simple and intuitive: If we switch on the 
coupling to the environment at a time when the system 
has already tunneled to the pair related island on the 
right (t = 56t), the environment will simply stop the 
state from tunnelling back to the initial island. Figure 
12 (right) shows the evolution of the Wigner function in 
this case. Contrariwise, if we turn on the coupling to 
the environment when the state is half way through the 
tunnelling (in an intermediate, delocalized, state at time 
to = 28t), decoherence yields an asymptotic state which 
has approximately half of the probability in each side of 
the well (this is seen on the left column of Figure 12). 
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FIG. 12. Evolution of the Gaussian state initially localized 
at the leftmost regular island (see Section 2A) when the en- 
vironment {D = 0.01) is connected at to = 28r (on the left) 
or at to = 56t (on the right). 
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VI. CONCLUSION 
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Quantum open systems become entangled with the en- 
vironments as a conscqiicncc of their interaction. For this 
reason, information initially stored in the state of the 
system irreversibly leaks into the environment. The re- 
sults we presented here support the point of view stating 
[14] that for classically chaotic systems the rate at which 
information flows from the system to the environment 
(the rate of von Neuman entropy production) is indepen- 
dent of the strength of the coupling to the environment 
(provided the coupling strength is above some threshold) . 
Thus, our results show that for chaotic quantum systems, 
the classical limit enforced by environment induced deco- 
herence is quite different from the one corresponding to 
regular systems. Thus, we showed that this limit exhibits 
an unavoidable source of unpredictability, being the rate 
at which information is lost into the environment entirely 
fixed by the chaotic nature of the Hamiltonian of the 
system. To the contrary, for regular systems the entropy 
production rate is proportional to the strength of the cou- 
pling to the environment. Therefore, the existence of this 
phase of coupling independent entropy production could 
be used as a diagnostic for quantum chaos. Our results 
also confirm previous estimations [14,17] for the transi- 
tion time tc between the two regimes: the one where the 
entropy production rate is diffusion-dominated and the 
one set by the chaotic dynamics, that is, a linear depen- 
dence of tc on the initial entropy (initial spread) and a 
logarithmic dependence on the diffusion constant. 

From our results it is possible to develop an intuitive 
picture for the reason why entropy production rate be- 
comes dominated by the chaotic dynamic. As we de- 
scribed in the paper, there are two interrelated processes 
contributing to the growth of entropy. First, the destruc- 
tion of the interference fringes which are dynamically 
produced in phase space by the streaching and folding 
of the chaotic evolution. The entropy production rate 
associated with this process is obviously determined by 
the slowest of the two timescales corresponding to the 
proceses of creation and anihillation of fringes. There- 
fore, when decoherence is effective, and the fringes diss- 
appear in a short timescale, the entropy production rate 
is of dynamical origin (and independent of the coupling 
strength between the system and the environment). On 
the other hand the spread of the possitive peaks of the 
Wigner function also contributes to the entropy growth. 
In this case, the rate becomes independent of the coupling 
strength once the phase space distribution approaches a 
critical width. Both of these processes are present in a 
general case (only it is possible to study them separately 
in idealized cases such as in the baker's map [46]). The 
results of this paper confirm this intuitive view, which 
can be made more precise when formulated in terms of 
equation (8). 

This work was partially supported by Ubacyt (TW23), 
Anpcyt, Conicet and Fundacion Antorchas. JPP thanks 
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